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Abstract 

Nonfactorizable effects in hadronic charmless B — > PP, VP decays can be parametrized in terms 
of the effective number of colors N^^ in the effective parameters a^^ that are linear combinations 
of Wilson coefficients. It is shown that N^^(V + A) in the penguin amplitudes induced by the 
{V — A){V + A) four-quark operators is different from N^^{V — A) in the decay amplitudes arising 
from the (V — A){V — A) operators. Central values of the branching ratios for uvr^ and 

B vrvr decays favor N^^(V — A) ~ 2, in accordance with the nonfactorizable effect observed 
in B ^ D^*'^tt(p). Measurements of the interference effects in B~ 7r~ (p~)'ir^{p^) decays will 
provide a more decisive test on the parameter N^^(V — A). However, N^^{V + A) ~ 2 is ruled 
out by 5=^ (j)K^. We find that the current bound on B^ 4>K^ imphes Nf{V + A) > 4.3, 
which is subject to the corrections from W^-annihilation and space-like penguin effects. With 
Nf{V -A) ^2 we show that the branching ratio of B ^ r]'K is enhanced considerably at small 
values of 1/N^^{V + A) so that it is compatible with experiment. In particular, the measurement of 

_^ rj' is now well explained without resorting to any new mechanism or new physics beyond 
the Standard Model. It is crucial to measure the charged and neutral decay modes of ^ (pK and 
B — > (j)K* to test the generalized factorization hypothesis. Finally, we point out that it is difficult 
to understand the observed large branching ratio of B^ — > ujK^ within the present framework. 
Inelastic final-state interactions may alleviate the difficulty with this decay mode. 
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I. INTRODUCTION 



To describe the hadronic weak decays of mesons, the mesonic matrix elments are cus- 
tomarily evaluated under the factorization hypothesis so that they are factorized into the 
product of two matrix elements of single currents, governed by decay constants and form 
factors. In the naive factorization approach, the relevant Wilson coefficient functions for 
color-allowed external W^-emission (or so-called "class-I" ) and color-suppressed (class-II) in- 
ternal 1^-emission amplitudes are given by ai = Ci + C2/Nc^ 02 = C2 + Ci/Nc, respectively, 
with Nc the number of colors. In spite of its tremendous simplicity, naive factorization en- 
counters two major difficulties. First, it never works for the decay rate of class-II decay 
modes, though it usually operates for class-I transition. For example, the predicted decay 
rate of the color-suppressed decay — > K^tt^ in the naive approach is too small when 
compared with experiment (for a review, see [1]). Second, the hadronic matrix element un- 
der factorization is renormalization scale n independent as the vector or axial-vector current 
is partially conserved. Consequently, the amplitude Cj(|u)(0)fact is not truly physical as the 
scale dependence of Wilson coefficients does not get compensation from the matrix elements. 
The first difficulty indicates that it is inevitable and mandatory to take into account non- 
factorizable contributions, especially for class-II decays, to render the color suppression of 
internal W emission ineffective. The second difficulty also should not occur since the matrix 
elements of four-quark operators ought to be evaluated in the same renormalization scheme 
as that for Wilson coefficients and renormalized at the same scale /x. 

Because there is only one single form factor (or Lorentz scalar) involved in the class-I or 
class II decay amplitude of B [D] PP, PV decays (P: pseudoscalar meson, V: vector 
meson), the effects of nonfactorization can be lumped into the effective parameters Oi and 
a2 [2]: * 



af = ci + C2( — + , af = C2 + ci (^— + X2) , (1-1) 



where Ci^2 are the Wilson coefficients of the spectator 4-quark operators, and nonfactorizable 
contributions are characterized by the parameters Xi and Xi- Taking the decay B~ D^7r~ 
as an example, we have [4-6] 

Xl — ^8 H ^1 ) X2 — ^8 H ^1 ) U-^J 



C2 Cl 



where 



(BD,.) ^ {D'n-\{du),_,{cb),_,\B-),,f ^ (D°7r-|(J^),_,(cb),_Ji?-) _ 

{D^7r-\{du),_,{cb),_,\B-)f {7r-\{du),_jO){D^\{cb),_jB-) ' 



*As pointed out in [3], the general amplitude of B{D) VV deeay consists of three independent 
Lorentz scalars, corresponding to S-, P- and D-wave amplitudes. Consequently, it is in general not 
possible to define an effective ai or 02 unless nonfactorizable terms contribute in equal weight to 
all partial wave amplitudes. 
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(BD,.) ^ 1 {D'7r-\idX'^u),_^ic\%),_jB-) 

are nonfactorizable terms originated from color-singlet and color-octet currents, respectively, 
{Qi(12)v-a = ^i7At(l - 75)^2, and (giA"g2)y_^ = giA"7^(l - 75)^2- The subscript 'f and 'nf 
in Eq. (1.3) stand for factorizable and nonfactorizable contributions, respectively, and the 
superscript {BD,7r) in Eq. (1.2) means that the pion is factored out in the factorizable 
amplitude of 5 — > Dtt and likewise for the superscript {Bti^D). In the large- A'^c limit, 
£1 = 0{1/N^) and = 0{1/N^ [6]. Therefore, the nonfactorizable term x in the Nc ^ 00 
limit is dominated by color octet-octet operators. Since |ci/c2| ^ 1, it is evident from Eq. (1) 
that even a small amount of nonfactorizable contributions will have a significant effect on 
the color-suppressed class-II amplitude. If xi,2 are universal (i.e. process independent) in 
charm or bottom decays, then we still have a new factorization scheme in which the de- 
cay amplitude is expressed in terms of factorizable contributions multiplied by the universal 
effective parameters a^^. (For B — > VV decays, new factorization implies that nonfactor- 
izable terms contribute in equal weight to all partial wave amplitudes so that can be 
defined.) The first systematical study of nonleptonic weak decays of heavy mesons within 
the framework of the generalized factorization was carried out by Bauer, Stech, and Wirbel 
[7]. Phenomenological analyses of two-body decay data of D and B mesons indicate that 
while the generalized factorization hypothesis in general works reasonably well, the effective 
parameters a'j;^ do show some variation from channel to channel, especially for the weak 
decays of charmed mesons [2,5,8]. An eminent feature emerged from the data analysis is 
that af^ is negative in charm decay, whereas it becomes positive in the two-body decays of 
the B meson [2,9,6]: 

af{D Kn) ~ -0.50 , af{B Dn) ~ 0.26 . (1.4) 

It should be stressed that since the magnitude of ai,2 depends on the model results for form 
factors, the above values of 02 should be considered as representative ones. The sign of af^ is 
fixed by the observed destructive interference in — > 'K^n'^ and constructive interference 
in B- D°n-. Eq. (1.4) then leads to 

X2il^ ~ rric, D — > Kn) ~ —0.36 , X2(a* ~ B — >• Dn) ~ 0.11 . (1.5) 

In general the determination of X2 is easier and more reliable than Xi- The observation 
|X2(-B)| -C |x2(-D)| is consistent with the intuitive picture that soft gluon effects become 
stronger when the final-state particles move slower, allowing more time for significant final- 
state interactions after hadronization [2]. 

Phenomenologically, it is often to treat the number of colors as a free parameter and 
fit it to the data. Theoretically, this amounts to defining an effective number of colors A^^*^, 
called l/i in [7], by 

l/ATf = (l/Arj+x. (1.6) 

It is clear from Eq. (1.5) that 
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(1.3) 



Nf{D Kti) » 3, Nf{B Dti) ^ 2. (1.7) 

Consequently, the empirical rule of discarding subleading l/N^ terms formulated in the 
large-A^c approach [10] is justified for exclusive charm decay; the dynamical origin of the 
1/Nc expansion comes from the fact that the Ficrz l/N^ terms are largely compensated by 
nonfactorizable effects in charm decay. Since the large-A^^c approach implies a'^ ~ C2 and 
since is observed to be positive in B~ — >• D^{t:~ , p~) decays, one may wonder why is the 
1/Nc expansion no longer applicable to the B meson ? Contrary to the common belief, a 
careful study shows this is not the case. As pointed out in [6], the large- A^c color counting 
rule for the Wilson coefficient C2(/x) is different at ~ and /i ~ rric due to the presence 
of the large logarithm at ~ rric. More specifically, C2{mb) — 0{1/N^ and C2{m^ — 0{1). 
RecaUing that Ci = (9(1), it follows that in the large- Ac limit [6]: 

f C2(mc) -I- 0(1/Ac) for the D meson, 

= |c2(m6) + ci(m6) (^ + £8(mb)) +0(1/A^c') for the S meson. ^^"^^ 

Therefore, a priori the 1/A^c expansion does not demand a negative for bottom decay ! 
and Nf^{B — > Dtt) ~ 2 is not in confiict with the large- Ac approach ! It should be remarked 
that although X2 is positive in two-body decays of the B meson, some theoretical argument 
suggests that it may become negative in high multiplicity decay modes [6]. 

Thus far the nonfactorizable effect is discussed at the purely phenomenological level. It 
is thus important to have a theoretical estimate of Xi even approximately. Unfortunately, 
all existing theoretical calculations based on the QCD sum rule [11], though confirm the 
cancellation between the 1/Ac Ficrz terms and nonfactorizable soft gluon effects [12], tend 
to predict a negative x in 5° — > D^7t~ , D^tt^ and B — > J/ipK decays. This tantalizing issue 
should be clarified and resolved in the near future. It is interesting to remark that, relying 
on a different approach, namely, the three-scale PQCD factorization theorem, to tackle the 
nonfactorizable effect, one of us and Li [13] are able to explain the sign change of X2 from 
bottom to charm decays. 

For B meson decay, the effective parameters af^2 have been determined so far only for 
B D{n,p) and B J/ipK where nonfactorizable effects amount to having A'c^ ~ 2. 
Recently, several exclusive charmless rare B decay modes have been reported for the first 
time by CLEO [14-19] and many of them are dominated by the penguin mechanism. It is 
thus important to know (i) does the constructive interference of tree amplitudes persist in 
class-Ill charmless B decay ? (class-Ill transitions receive contributions from both external 
and internal W emissions) , and (ii) is A^c''^ same in spectator and penguin amplitudes ? In 
the literature it is customary to assume that A^^ behaves in the same way in the penguin and 
non-penguin amplitudes. The decay rate of the rare B decays is then studied as a function of 
1/A'c''^. In the present paper, we shall sec that, theoretically and experimentally, N^^ (V + A) 
in the penguin amplitude induced by the {V — A){V + A) quark operators is different from 
Nf^{V — A) in the tree or penguin amplitude induced by the {V — A){y — A) operators. ^ 



^In [20] we have assumed that N^^{V — A) ^ Nf^{V -|- A) 2. The present paper is an improved 
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We find that N^^{V + A) in penguin- dominated charmless B decays is clearly larger than 
Nf^[V — A) extracted from spectator-dominated processes. Therefore, the nonfactorizable 
effect in tree and penguin diagrams behaves in a different manner. This observation is the 
key element for understanding the CLEO measurement of B"^ — > -q'K"^ and B^ rj'K°. By 
treating Nf{V - A) and Nf{V + A) differently, the data of S ^ rj'K can be explained in 
the framework of the Standard Model without resorting to new mechanisms or new physics 
beyond the Standard Model. 

This paper is organized as follows, in Sec. II we sketch the starting point of the effective 
Hamiltonian and emphasize that vertex and penguin corrections to the four-quark operators 
should be combined together with the Wilson coefficients to render the resulting physical 
amplitude independent of the choice of the renormalization scheme and scale. Then we 
extract the information of Nf^{V — A) from spectator dominated charmless B decays B"^ 
Lon^ and S — > tttt in Sec. Ill and Nf^{V + A) from the penguin dominated process B — > (f)K 
in Sec. IV. In Sec. V we demonstrate that the measurement of B ^ rj'K also favors a 
different treatment of N^^iV — A) and N^^{V + A). In Sec. VI we point out a difficulty 
with B"^ — * ujK"^ within the present framework. Conclusions and discussions are presented 
in Sec. VII. 



II. CALCULATIONAL FRAMEWORK 



We briefly sketch in this section the calculational framework. The relevant effective 
AB — 1 weak Hamiltonian is 
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Kff(As = 1) = ^ [K6-K;(cior + C2O2") + Vcbv:j{c,ot + c^o^^ - Vt,v:^ ^ + h.c, (2.1) 

where q — d,s, and 



03(5) = (#)y-AlIW)v-A(v+A)> 



C7(9) = ^{qb)y_^J2^l'W9')v+Aiv-A), 



01 



04(6) = {qabp)y_^ Y.Wa) 



v-a{v+a)^ 



(2.2) 



08(10) = 7;i(labi3)y_^J2eg,{qpqJ 



V+A 



(v-a)' 



with O^-Oq being the QCD penguin operators and O7-O10 the electroweak penguin operators. 
As noted in passing, in order to ensure the renormalization-scale and -scheme independence 
for the physical amplitude, the matrix elements of 4-quark operators have to be evaluated 
in the same renormalization scheme as that for Wilson coefficients and renormalized at the 
same scale /i. 

In full theory, the leading QCD correction to the weak transition is of the form 
as ln(M^ fo^ massless quarks, where p is the off-shell momentum of external quark lines 



version of [20]. 
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and its magnitude —p^ depends on the system under consideration. For example, —p^ ~ ml 
in the energetic two-body charmless B decays. The merit of the effective Hamiltonian ap- 
proach is that one can choose a renormalization scale so that the leading logarithmic 
correction ln(M^/ — p^) = ln(M^//x^) -|- ln(//^/ — p^) is decomposed in such a way that 
the large logarithmic term ln(M|^///^) is lumped into the Wilson coefficient function c(//) 
and summed over to all orders in ag using the renormalization group equation, while the 
logarithmic correction ln(/i^/— p^) to the matrix element (0(/i)) is small (for a review, see 
[26]). Since 0(/i) is the four-quark operator renormalized at the scale /x, its hadronic matrix 
element is related to the tree level one via 



{0{^X)) = ^(//)(0)tree, 



(2.3) 



with 



(2.4) 



for current-current operators, where we have included the non-logarithmic constant contri- 
bution c since the logarithmic contribution ln(//^/— p^) is small when /x^ ~ —p^ and hence the 
momentum-independent constant term cannot be neglected. It follows that schematically 



tree 



{O) 



tree- 



(2.5) 



To the next-to-leading order (NLO), c(/x) depends on the renormalization scheme chosen, so 
does the constant c in g{ij). However, the effective Wilson coefficient c^^ is independent of 
the choice of the renormalization scheme and scale. It should be stressed that, except for the 
lattice QCD, model calculations of the hadronic matrix elements are actually performed for 
(O)tree rather than for (0(/x)). (Quark model calculation of (O)tree, for example, may involve 
an implicit low energy scale, but it has nothing to do with the renormalization scale /x.) For 
example, in the factorization approximation, the matrix element (O)fact is scale independent 
and hence it cannot be identified with (0(//)). Therefore, it is important to evaluate g{^x)i 
the perturbative corrections to the four-quark operators at the scale 

As emphasized above, before applying factorization or carrying out any model calculation 
of hadronic matrix elements, it is necessary to incorporate QCD and electroweak corrections 
to the operators: 



so that Ci{ii){Oi{ii)) = cf{Oi) tree, where 



I+^^^(/^) + i^^e(/x) 



tree) 



(2.6) 



.eff 



(2.7) 



Then the factorization approximation is applied to the hadronic matrix elements of the 
operator O at tree level. Perturbative QCD and electroweak corrections to the matrices 
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rhs and rhe from vertex diagrams and penguin diagrams have been calculated in [21-24]. 
It should be remarked that although the penguin coefficients C3 — Cio are governed by the 
penguin diagrams with t quark exchange, the effective Wilson coefficients do incorporate the 
effects of the penguin diagrams with internal u and c quarks induced by the current-current 
operator Oi. For example [24], 



Stt 



cM + ---, (2.8) 



where Aj = VibV*^ {q = d,s), G{mq, k, n) = — G{'mq, k,n), k is a constant depending on 
the renormalization scheme, k is the gluon's virtual momentum, and 

G{m, k, 11) = -4 dxx{l - x) In - - ^2.9) 

For b ^ s transitions, |Au| <^ |At|, Ac ~ —At, and hence 

cf = C6(/i) + ^G(m„ k, /x)ci(/x) + • • • . (2.10) 

The importance of the so-called "charming" penguins for 6 — > s transition was emphasized 
recently (and probably over-emphasized) in [25]. 

Using the next-to-leading order AB = 1 Wilson coefficients obtained in the 't Hooft- 
Veltman (HV) scheme and the naive dimension regularization (NDR) scheme at // = m5(m{,), 
Ag^ = 225 MeV and mt = 170 GeV in Table 22 of [26], we obtain the effective 
renormahzation-scheme and -scale independent Wilson coefficients cf^ at k'^ — ml/ 2: ^ 

cf = 1.149, cf = -0.325, 

cf = 0.0211 + i0.0045, cf = -0.0450 - i0.0136, 

cf = 0.0134 + i0.0045, cf = -0.0560 - i0.0136, 

cf = -(0.0276 + i0.0369)Q;, cf = 0.054 a, 

cf = -(1.318 + i0.0369)Q;, cf = 0.263Q;. (2.11) 

Two important remarks are in order. First of all, c^^ are surprisingly very close to the leading 
order Wilson coefficients: c\'^ — 1.144 and Cs*^ = —0.308 &t /i — mb{mb) [26], recaUing that 
^NDR ^ _o.i85 and cf^ = -0.228 at NLO [26] deviate substantially from the leading order 
values. This means that (Oi,2(a*)) ~ (Oi^2)tree- Hence, it explains why the conventional way 
of applying the Wilson coefficients at leading order and evaluating the matrix elements of 
current-current operators at tree level is "accidentally" justified provided that /x^ ~ —p^. 
Second, comparing (2.11) with the leading-order penguin coefficients [26] 



■tWe use the complete expressions of rhs{iJi) given in [24] and rhein) in [22] to evaluate cf . Note 
that while our cf^ are consistent with the numerical results given in [24,27], our values for cf 
are different from that shown in [27]. 
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c^" = 0.014, = -0.030, df = 0.009, c^" = -0.038 (2.12) 

ji — mb{mb), we see that Re(c|^g) ^ ^Cs-eip')- This imphes that, contrary to the case 
of current-current operators, penguin corrections to the current-current operators give im- 
portant contributions to the QCD penguin operators. This means that the decay rates of 
charmless B decay modes dominated by penguin diagrams will be too small by a factor of 
~ (1-5)^ = 2.3 if only leading-order penguin coefficients are employed for calculation. 

We shall see later that running quark masses appear in the matrix elements of {S — 
P) {S + P) penguin operators through the use of equations of motion. The running quark 
mass should be applied at the scale // ~ m;, because the energy release in the energetic 
two-body charmless decays of the B meson is of order nib. Explicitly, we use [28] 

mu{mb) = 3.2 MeV, md{mb) = 6.4 MeV, ms{mb) = 105 MeV, 

mdmb) = 0.95 GeV, mb{mb) = 4.34 GeV, (2.13) 

in ensuing calculation, where we have applied = 150 MeV at /i — 1 GeV. 

It is convenient to parametrize the quark mixing matrix in terms of the Wolfenstein 
parameters: A, A, p and r] [29], where A — 0.804 and A = 0.22. In the present paper we 
employ two representative values for p and 77: (i) p = 0.16, r) = 0.34, and (ii) p = —0.12, rj = 
0.35. Both of them satisfy the constraint \/ p'^ + 1]'^ = 0.37. A recent analysis of all available 
experimental constraints imposed on the Wolfenstein parameters yields [30] 

p= 0.156 ±0.090, = 0.328 ± 0.054, (2.14) 

where p = p{l — ^) and f] = r]{l — and it implies that the negative p region is excluded 
at 93% C.L.. 



III. NONFACTORIZABLE EFFECTS IN SPECTATOR AMPLITUDES 

The combinations of the effective Wilson coefficients a2i = d§ ]^c|f_i, a2i-i = c|f_i -|- 
W'^'^ ~ 1) ■ ■ ■ ) 5) appear in the decay amplitudes. As discussed in the Introduction, 
nonfactorizable effects in the decay amplitudes of S — > PP, VP can be absorbed into the 
parameters af^ . This amounts to replacing A^^ in Oj by {Nf^)i. (It must be emphasized 
that the factor of Nc appearing in any place other than should not be replaced by Nf^.) 
Explicitly, 

eflf _ cflf , 1 „cflF cff _ cflf , ^ pCff fQ D 

It is customary to assume in the literature that (A^^)i {Nf^)2 ■ ■ ■ {NI^)iq so that the 
subscript i can be dropped. A closer investigation shows that this is not the case. Consider 
an operator of the form O — F^f Q3^'q2 which arises from the Fierz transformation of a 
singlet-singlet operator with F and F' being some combinations of Dirac matrices. Applying 
the identity 
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O = ^?ir?2 qs^'q^ + ^giAT?2 ?3AT'g4, (3.2) 
to the matrix element of M — > P1P2 leads to (assuming the quark content 51^2 for Pi) 

(P1P2IOIM) = ^{P,\q,rq2\0){P,\qsr'q,\M) + ^{PiP2\qirq2q3r'q^\M)^f 

+ ^(PiP2|giATg2g3A"r'g4|M). (3.3) 

The nonfactorizable effects due to octet-octet and singlet-singlet operators are characterized 
by the parameters £8 and ei, respectively, as shown in Eq. (1.3): 

_ {P,P2\qiTq2q3T'q,\M)^f {PiP2\qiTq2 qs^'q^lM) f 

^1 ~ , „ ^ , ^ — — r- — — — — X 



{PiP2\qirq2q3r'q,\M)f (P1P2I (giga),^^ (g3g4)v_ JM)/ 
^ 1 (PiP2|giA"rg2g3A"rg4|M) {P,P2\q,Tq2 gsT' q4\M) j 

2 {P,P2\q,Tq2q3T'q,\M)f {P,P2\{qiq2)y.AM4)y.jM)/ 

However, the Fierz transformation of the {V — A){V + A) operators 05,6,7,8 is quite different 
from that of {V — A){V — A) operators Oi,2,3,4 and 09,10; that is, 

{V -A){V + A)^ -2{S -P){S + P), 

{V-A){V-A)^{V-A){V-A). (3.5) 

Therefore, F and T' arc the combinations of the Dirac matrices 1 and 75 for the Fierz 
transformation of {V — A){V + A) operators, and the combinations of 7^ and 7^75 for 
{V — A){V — A) operators. As a result, nonfactorizable effects in the matrix elements of 
[V — A){V + A) operators are a priori different from that of {V — A){y — A) operators, i.e. 
X{V + A)^ x{V - A). Since ^l/N^ + x [cf. Eq. (1.6)], theoretically it is expected 

that 



Nf{v -A)^ {n:^)_^ ^ (ivf )^ ^ (n:^),^ ^ (n:^)^ « (ivf )^ ^ (ivf , 

Nf{V + A) = (Nf)^ ^ (Nf)^ ^ (Nf)^ ^ (Nf)^ , (3.6) 

and N^^(y+A) 7^ N^^{V—A) in general. In principle, A''^*^ can vary from channel to channel, 
as in the case of charm decay. However, in the energetic two-body B decays, A*"^''^ is expected 
to be process insensitive as supported by data [6]. As stressed in the Introduction, if A^^ is 
process independent, then we have a generalized factorization. Contrary to the naive one, 
the improved factorization does incorporate nonfactorizable effects in a process independent 
form. For example, xi = X2 = — | in the large- Ac approximation of factorization. 

The unknown parameters {Nf^)i in charmless B decays in principle can be determined 
if the decay rates are measured for a handful of decay modes with sufficient accuracy. Due 
to the limited data and limited significance available at present we shall use (3.6) and the 
experimental result for N^^{B — > Dn) as a guidance to determine {N^^)^. To begin with, we 
focus in this section the decay modes dominated by the spectator diagrams induced by the 
current-current operators Oi and O2. In particular, we would like to study these modes which 
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arc sensitive to the interference between external and internal 1^-emission amplitudes. The 
fact that N^^ < 3.5 {Nf^ > 3.5) implies a positive (negative) a"^ and hence a constructive 
(destructive) interference will enable us to differentiate between them. Good examples are 
the class-Ill modes: B"^ — > un^, n^n"^, 7771=*=, n^p^, • • •, etc. 

We first consider the decay B" — > u;7r~. Under the generahzed factorization, its decay 
amplitude is given by 



td 



04 + aio - 2(a6 + as) 



ml 



{nib + mu){mu + rrid) J 



+ ^ (4a3 + 2a4 + 4a5 + 07 + Og - aio) X^^""'"^^ 



+ 2 ( 04 + oio - 2(a6 + 08 ) 



(mb + m„)(m„ + m^) / 



(3.7) 



where we have dropped the superscript "eff" for convenience, and the notation for 
example, denotes the factorization amplitude with the tt meson being factored out: 



= (o;|(««),_jO)(7r-|(j6),_Ji?-), 
^ (7r-(^|(J«),_jO)(0|(^6),_j5-) 



(3.8) 



Note that in the penguin amphtude, the term X^^''^'^^ arises from the space-like penguin 
diagram. Using the following parametrization for decay constants and form factors: ^ 



(0|A^|P(g)) = ifpq^,, (0|V;|\/(p,£)) = /ymv£^. 



(p'(p')|v;.|p(p)) = (pm+p;.- ^2 



gjFi(g^) + Fo(g^) 



2 2 
2^ TTlp TTl pi 



{V{p',e)\V,\P{p)) 



mp + mv 



{V{p\e)\A^\P{p))^i {mp + mv)e^A,{q^) - 



e ■ p 



mp + mv 



e • p 



2mv^qM^i.r)-Mr 



(3.9) 



where q = p-p', Fi(0) = Fo(0), ^3(0) = Ao(0), and 



§Once the one-body matrix elements are defined, one can apply heavy quark symmetry to the 
two-body matrix elements for heavy-to-heavy transition to show that all the form factors defined 
in (3.9) are positively defined at > and that the relative signs between two-body and one-body 
matrix elements are fixed. In this way, we find that the vector form factor V^(g^) defined by Bauer, 
Stech and Wirbel [31] has a sign opposite to ours. Note that our convention is eoi23 = 1- 
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A{q ) = -^;^M<1 ) - (3-10) 

we obtain 

XiBn,-) ^ -V2Um^Ff^ml)e ■ ps. (3.11) 

For the dependence of form factors in the region where (f is not too large, we shall 
use the pole dominance ansatz, namely, 

= n 2^n , (3.12) 

(1 - q^/ml) 

where is the pole mass given in [7]. A direct calculation of B ^ P and B ^ V form 
factors at time-like momentum transfer is available in the relativistic light-front quark model 
[32] with the results that the q"^ dependence of the form factors A^, Fi is a dipole behavior 
(i.e. n — 2), while Fq exhibits a monopole dependence {n — 1). The decay rate is then given 

by 



V / o 2 



A{B- TT-cu) 



£ -PB 



(3.13) 



where Pc is the cm. momentum 



J[ml - (^o; + m^y][ml - (m^ - m^Y] 
2^, ■ 

Since 

Vu,V::, = AX%p-t7i), Vc,V:,^-AX\ Vt,V;, = AX%l-p + t7i), (3.15) 

in terms of the Wolfenstein parametrization [29], are of the same order of magnitude, it is 
clear that B~ — > U7r~ is dominated by external and internal W emissions and that penguin 
contributions are suppressed by the smallness of penguin coefficients. Neglecting the W- 
annihilation contribution denoted by X^^''^'^\ and using = 132 MeV, = 195 MeV for 
decay constants, A^'^(O) = 0.28/^, Ff"(0) = 0.33 for form factors [31], and t{B^) = 
(1.67 ± 0.04) ps [33] for the charged B hfetime, the branching ratio of B^ n^uj averaged 
over CP-conjugate modes is shown in Fig. 1 where we have set N^^{V + A) = N^^{V — ^4) = 
Nf^ and plotted the branching ratio as a function of We see that the branching 

ratio is sensitive to l/Nf^ and has the lowest value of order 2 x 10~^ at Nf^ — oo and then 
increases with 1/Nf^. Since experimentally [16] ** 



**The significance of B ujtt is reduced in the recent CLEO analysis and only an upper limit 
is quoted [34,19]: B{B^ tt^uj) < 2.3 x lO'^. Since B{B^ K^ui) = (1.5+g:^ ± 0.2) x IQ-^ and 
B{B^ h^oj) = {2.5to^ ± 0.3) x 10"^ with h = Tr, K, the central value of B{B^ tt^lo) remains 
about the same as (3.16). 
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B{B^ UTT^) = (l.llSj ± 0.2) X 10-^ (3.16) 

it is evident that l/Nf^ > 0.35 is preferred by the data. Because this decay is dominated by 
tree amphtudes, this in turn imphes that 

Nf{V -A) < 2.9 from n^u. (3.17) 

With the value of Nf^{V — A) being fixed to be 2, the branching ratio of B"^ — > tt^o; is 
plotted in Fig. 2 as a function of Nf^{V + A). We see that for positive p, which is preferred 
by the current analysis [30], the branching ratio is of order (0.9 — 1.0) x 10~^, very close to 
the central value of the measured one. 

The fact that N^^{V—A) < 2.9 in charmless two-body decays of the B meson is consistent 
with the nonfactorizable term extracted from B {D,D*)7r, Dp decays, namely x ~ 0-10 
or N^^{B — >■ Dtt) ~ 2. Since the energy release in the energetic two-body decays B lvtt, 
B — > Dtt is of the same order of magnitude, it is thus expected that N^^{V — A)\b^uj77 ~ 2. 

The main uncertainty of the above analysis is the negligence of the space-like penguins 
and VF-annihilations. It is common to argue that 1^- annihilation is negligible due to helicity 
suppression, corresponding to form factor suppression at large momentum transfer, = m% 
(for a recent study, see [35]). However, we see from Eq. (3.7) that the space-like penguin 
contribution gains a large enhancement by a factor of m^/[m5(m„ + 171^)] ~ 670. Therefore, 
there is no good reason to ignore the space-like penguin effect [36] that has been largely 
overlooked in the literature. Unfortunately, we do not have a reliable method for estimating 
1^- annihilation and hence space-like penguins. 

We next come to the decay B~ — > 7r~7r° which is quite clean and unique in the sense that 
this is the only two-body charmless B decay mode that does not receive any contributions 
from the QCD penguin operators. Under the generalized factorization approximation, 

A{B- tt-tt") = ^V^,V:,{a^ + a,)iU{ml - ml)Ft\ml), (3.18) 

with Fq'"° = Fq'"^ I where we have neglected the very small electroweak penguin contri- 
butions. The decay rate is 

r(B- ^ 7r-7r°) = \MB~ ^ 7r-7r°)P. (3.19) 

871171% 

Just like the decay B~ 7r~cj, the branching ratio of B~ — > 7r~7r° also increases with 
as shown in Fig. 3. The CLEO measurement is [17] 

B{B^ 7r±7r°) = (o.9^°j) x 10"^ < 2.0 x 10"^ (3.20) 

However, the errors are so large that it is meaningless to put a sensible constraint on Nf^iV — 
A). Nevertheless, we see that in the range < 1/Nf < 0.5 [24], Nf{V - A) ^ 2 is favored. 

In analogue to the decays B D^*^n{p), the interference effect of spectator amplitudes 
in class-Ill charmless B decay can be tested by measuring the ratios: 
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Since penguin contributions are very small as we have checked numerically, to a good ap- 
proximation we have 




U A^%ml) aX 

^3-^(l + -^^-J . (3.22) 

Evidently, the ratios Ri are greater (less) than unity when the interference is constructive 
(destructive). Numerically we find 

z?_/l-74, z?_/l-40' / 2.50 for iVf = 2, 



0.58, ^ 10.80, [0.26 for A^f = oo, 

where use of t{B^) ^ (1.57 ± 0.04) ps [33], fp = 216 MeV, ^^''(0) = 0.28 [32] has been 
made. Hence, a measurement of (in particular R^), which has the advantage of being 
independent of the parameters p and rj, will constitute a very useful test on the effective 
number of colors N^^{V — A). The present experimental information on S — > tt'^tt" is [17] 

B(B° TT^TT^) = (0.7 ± 0.4) X 10"^ < 1.5 x 10"^ (3.24) 

As far as the experimental central value of Ri is concerned, it appears that l/Nf^ ~ 0.5 is 
more favored than any other small values of . 

In short, using the central values of the branching ratios for class-Ill decay modes: B 
TTUJ, B irn, wc find that within the range < < 0.5, N^^{V — A) ~ 2 is certainly 

more preferred. Measurements of class-Ill decays are urgently needed in order to pin down 
the nonfactorizable effect in tree amplitudes. In particular, measurements of the interference 
effects in charged B decays B~ — > 7r~(p~)7r'^(p°) will be very helpful in determining Nf^{V — 
A). 



IV. NONFACTORIZABLE EFFECTS IN PENGUIN AMPLITUDES 

In Sec. Ill we have shown that for spectator-diagram amplitudes Nf{V - A) ~ 2 is 
preferred, as expected. However, the nonfactorizable effect in the penguin amplitude is not 
necessarily the same as that in the tree amplitude since the chiral structure and the Fierz 
transformation of the (V — A)[V + A) 4-quark operators 05,6,7,8 are different from that of 
(V - A){V - A) operators; that is, Nf{V + A) is a priori not the same as Nf{V - A). 
By studying the penguin-dominated decays B — > (f)K and B — > (f)K*, we shall see that 
Nf^{V -I- A) ~ 2 is ruled out by the current bound on B~ — > K~(f). 



13 



The decay B 
A{B- 



K (f) receives contributions from VT-annihilation and penguin diagrams: 



as + 04 + as - ^(ay + 09 + Oio)) X^-^^'"^) 



+ 04 + oio - 2(a6 + as) 



ml 



(mb + mu){ms + 



j^(B,K4,) 



(4.1) 



where 



\B- 



(4.2) 



Neglecting 1^-annihilation and space-hke penguin diagrams and using = 237 MeV, 
^(0) = 0.34 [32], we plot in Fig. 4 the branching ratio of B^ (t)K^ against 1/A^^^*^ for 
two different cases: the dotted curve for the free parameter Nl^(y-\-A) — Nl^(y — A) — Nf^ 
and the solid curve with Nf^iV — A) being fixed at the value of 2. In either case, it is clear 
that N^^{V + A) = 2 is evidently excluded from the present CLEO upper limit [34] 



B{B^ (j)K^) < 0.5 X 10 



(4.3) 



A similar observation was also made in [37]. The conclusion that N^^{V + A) ^ 2 will be 
further reinforced if the decay rate of B^ — > (j)K^ is enhanced by the space-like penguins. 
From Fig. 4 we also sec that l/Nf{V + A) < 0.23 or Nf{V + A) > 4.3. Note that 
this constraint is subject to the corrections from space-like penguin and VT-annihilation 
contributions. At any rate, it is safe to conclude that 



Nf{V + A)> Nf(V - A). 



eff/ 



(4.4) 



The branching ratio of B ^ <pK*, the average of (pK* and (pK*^ modes, is also measured 
recently by CLEO with the result [34] 



B{B (f)K*) = ^ [B{B^ (t>K*^) + B{B^ 0^*°)] = (l.ll^l ± O.2) x 



10" 



(4.5) 



As emphasized in the first footnote of Sec. 1, the effective parameters ai in general cannot be 
defined for the B — > VV decay as its amplitude involves more than one Lorcntz scalar. In 
the absence of information for the nonfactorizable contributions to various Lorentz scalars, 
we shall assume generalized factorization. Under this hypothesis, the decay amplitude of 
B — > (t)K* has a similar expression as that of S — > (t)K. Its decay rate is given by 



T{B- (j)K*-) 



Pc 



G 



X ( 03 04 05 - - (ay Og Oio) 



(a-6x)' + 2(l + cV) , (4.6) 



with 
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(4.7) 



where we have neglected contributions proportional to X^^'^*'^\ 

We calculate the decay rates using two different sets of values for form factors: 



Af^*(0) = 0.328, 



Af^ (0) = 0.331, 



y«^*(0) = 0.369 



(4.8) 



from [7] and 



^f^*(0) = 0.26, 



Af^*(0) = 0.23, 



T/S^*(0) = 0.32 



(4.9) 



from [32]. As for the dependence, light-front calculations indicate a dipole behavior for 
V{q'^), A2{q'^) and a monopole dependence for Ai(g^) [32]. The result is shown in Fig. 5. It 
is interesting to note that the branching ratios are very insensitive to the choice of the values 
for form factors, (4.8) or (4.9). We see that the allowed region is 0.7 > l/Nf{V + A) > 0.25 
or 4 ^ Nf^{V + A) ^ 1.4, bearing in mind that this constraint is subject to the corrections 
from annihilation terms. This seems to be in contradiction to the constraint Nf^{V + A) > 
4.3 derived from (j)K^. In fact, it is expected in the factorization approach that 

r{B ^ (j)K*) T{B ^ (j)K) when the W^-annihilation type of contributions is neglected. 
The current CLEO measurements (4.3) and (4.5) arc obviously not consistent with the 
prediction based on factorization. One possibility is that generalized factorization is not 
applicable to -B ^ VV . Therefore, the discrepancy between B{B — cfyK) and B{B <pK*) 
will measure the degree of deviation from the generalized factorization that has been applied 
io B ^ (f)K*. At any rate, in order to clarify this issue and to pin down the effective number 
of colors Nf^{V + A), we need measurements of S — > (f)K and B — > (f)K*, especially the 
neutral modes, with sufficient accuracy. 

Since CLEO has measured B^ — > n'K and B 7r~^K~ [17], we have also studied these 
two decay modes. We found that for a fixed N^^{V — A) = 2, the predicted branching ratios 
of i? TiK are in agreement with the CLEO measurement within errors for all values of 
1/N^^(y + A). Hence, no useful constraint on N^^iV + A) can be derived from B^ — > ir^K 
and S° ^ n+K-. 

Since Nf^{V + A) > Nf^{V — A), one may wonder if the leading 1/A^c expansion may 
happen to be applicable again to the matrix elements of {V — A){V + A) operators. We 
believe that Nf{V + A) = cxD is very unlikely for two reasons. First, it will predict a too 
small branching ratio of -B — (f)K* as shown in Fig. 5. Second, it implies a nonfactorizable 
term xb{V + A) ~ — |, as in the charm case. Since the energy release in the energetic 
two-body decays of the B meson is much larger than that in charm decay, it is thus expected 
that 



Because Nf{V + A) > 4.3, it is then plausible to assume that xsiV + A) ~ -x{B Dn) fa 
-(0.10 - 0.12). Hence, Nf{V + A) ^ (4.3 - 4.9). 



\X{D ^ Kn)\ > \xb{V + A)\ ~ \x{B ^ Dn)\. 



(4.10) 
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V. IMPLICATIONS ON CHARMLESS B DECAYS INTO rj' AND -q 



When the prehminary CLEO measurement of — > rj'K^ was reported last year [14] 

B{B^ Tj'K^) = (j.Stli ± I.O) X 10-^ (5.1) 

it has stimulated a great interest in the community since early theoretical estimates of the 
S± r]'K^ branching ratio [36,38,23] he in the range of (1 - 2) x 10"^^ Since then, 
many theoretical studies and speculation have surged, as evidenced by the recent literature 
[24,27,39-49] that offer various interpretations on the abnormally large branching ratios. It 
was soon realized [24,27] that the running strange quark mass at the scale ^ = 0{mh) and 
SU(3) breaking in the decay constants of the rjo and rj^ will provide a large enhancement 
to the decay rate of — > rj'K (for a review, see [50]). Unfortunately, as pointed out in 
[24], this enhancement is partially washed out by the anomaly contribution to the matrix 
element (77'|s75s|0), an effect overlooked previously. As a consequence, the branching ratio of 
B — > rj'K is of order (2 — 3) x 10~^ in the range < < 0.5. The discrepancy between 

theory and current measurements [18] 

B{B^ rj'K^) = {Q.btli ± 0.9) x 10-^ 

B{B^ ^ Vi^°) = (4.7^1^ ± 0.9) X 10-^ (5.2) 

seems to call for some new mechanisms unique to the rj' production or even some new physics 
beyond the Standard Model. 

All the previous analyses of B ^ rj'K in the literature are based on the assumption that 
{Nf^)i are the same for i = 1, 2, ■ ■ ■ , 10. In this Section we will show that the fact that 
N^^iy + A) and N^^{V — A) are not the same and that they are subject to the constraints 
(3.17) and (4.10) will lead to a significant enhancement for the decay rate of B ^ rj'K 
at small values of l/Nf^. Moreover, we shall see that the prediction of B{B — > rj'K) is 
compatible with experiment. Especially, the measurement of S° — > rj'K'^ is well explained, 
implying that no new mechanism in the Standard Model or new physics beyond the Standard 
Model is needed to account for the data. 

To begin with, we write down the factorizable amplitude 

A{B- ^ n'K-) = {a^X(^^''^^ + a^Xf ^'''') + a.X^^^^'^^) + V^V:,a,Xi^^^^'^ 



a, + aio + 2(a6 + a,)- —-^ -] 

^ [ms + m„) [mi, - m„) J 



ttThe prediction B{B^ r]'K^) = 3.6 x 10^^ given in [36] is too large by about a factor of 2 
because the normalization constant of the rj' wave function was not taken into account in the form 
factor Fq ' . This negligence was also erroneously made in some recent papers on 5 — > rj K. Note 
that all early calculations [36,38,23] did not take into account the anomaly contribution to the 
matrix element (77'|s75s|0) (see below). 



16 



+ 04 + aio + 2(a6 + ag) 



/ 111 
+ I as + a4 - a5 + -ar - -ag - -aio 



1 N 

+ (ae - Tragj 



1 _ 51 ) W(Bi^,'7') 



for B — > ?7'K , and 



2 ms{mb-ms) \ 



(5.3) 



+ ( 2a3 - 2a5 



4 



1 777^ 

a4 - 7^«io + (2a6 - a8)7 ■ r, 

2 {rris + md){mb - md) J 



(2as - 2a, - ^ar + ^ag) ^''''^ + (aa - a, - a^ + agjXf ^'''') 



+ a4 + aio + 2(a6 + ag) 



m, 



(rris - md)(mb + md) J 



>'K) 



' 111 

03 + 04 — as + -aj — -ag — 2^10 



1 N 

+ (ae - Tras) 



2 ms{mb-ms) \ f^, 
for — > r^'TC , where 

^ (X-|(,-^)^_jo)(r/'|(«6),_jS-) = (:^|(.-d),_jO)(r^'|(J6),_jS° 

^'^') ^ (77l(?-g),_JO)(X-|(.-6),_j5-) = (,7'|(5-g),_JO)(F°|(.-fe),_jS°) 
X^B,,'K) ^ (^'X-|(.-^),_J0)(0|(^6),_JS-), 



(5.4) 



(5.5) 



and use of the isospin relation xf"''''^ = has been made. For the amphtude of 

B~ r]'K^ , the terms proportional to X'^^'^'^^ and Xj^^^'^'^ with penguin coefficients are 
often missed or not considered in previous analyses. Note that the neutral mode b''' — rj'TC^ 
differs from the charged mode that it does not receive contributions from external l^-emission 
and VF-annihilation diagrams. Prom the relevant quark mixing angles 



VtbVZ = -AX' + ^Ail - 2p)\^ + ii^A\\ 



(5.6) 



it is clear that B — > r]'K decays are dominated by penguin diagrams. 
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The presence of the term 1 — (Z^"//,^/) in (5.3) and (5.4) is necessary and mandatory in 
order to ensure a correct chiral-limit behavior for the {S — P){S + P) matrix elements of the 
penguin operators 05,6,7,8- In the chiral limit m^, ma, rus 0, the ratio m\/{ms + my) — 
m^/ (m„ + md) remains finite-'--'-, but this is no longer the case for m^,/ms associated -with the 
matrix element {r)'\sj5s\0) since the 77' mass originates from the QCD anomaly and does not 
vanish in the chiral limit. As pointed out in [42,24], due to the presence of the anomaly in 
the equation of motion 

5'^(s7m75s) = 2m,si75S + -^0^.0"^", (5.7) 
it is erroneous to apply the relation 

2 

WMO)^-t^J^„ (5.8) 

as adopted previously in the literature, -where (0|g7^75g|r7') = if^iPn- Neglecting the u and 
d quark masses in the equations of motion leads to [51] 

(r7'|gGG|0) = /X' (5-9) 

and hence [42,24] 

2 

(r;'k-75^|0) = -i^ {r,, - f:;,) . (5.10) 



s 



It is easily seen that this matrix element has the correct chiral behavior. It should be stressed 
that in order to go the chiral-symmetry limit, one must consider both and ^ — 

together [52], where 9 is the rj — r]' mixing angle to be defined below. Since f^, ^ \f^i 
(see below) and the decay amplitude is dominated by {S — P){S + P) matrix elements, it 
is obvious that the decay rate of 5 — > 'q'K is reduced considerably by the presence of the 
anomaly term in (r;'|s75s|0). 

To determine the decay constant , we need to know the wave functions of the physical 
77' and 77 states which are related to that of the SU(3) singlet state 770 and octet state 7/8 by 

77' = 778 sin ^ + 770 cos 6*, 77 = 778 cos ^ — 770 sin ^, (5-11) 

with 9 —20°. When the rj — r]' mixing angle is —19.5°, the 77' and 77 wave functions have 
simple expressions [36]: 

\r]') — --j=\uu + dd + 2ss) , \r]) — -j=\uu + dd — ss) , (5-12) 

recalling that 



■'■■'■For the annihilation term, the chiral-Umit behavior of — r—^ x^^''^'^'^ is supposed to be 

taken care of by the form factors in X^^''^'^\ 
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\rjo) — —^\uu + dd + ss), \rjs) — —7=\uu + dd — 2ss). (5.13) 

At this specific mixing angle, f^, — |/*/ in the SU(3) hmit. Introducing the decay constants 
fs and /o by 

then and /*/ are related to fs and /o by 

P, = ^ sin ^ + A cos ^, = -2 A sin ^ + A cos ^. (5.15) 

Likewise, for the rj meson 

= ^ cos ^ - A sin P = -2 A cos ^ - A sin 0. (5.16) 

The factorizable amplitude denoted by Xj^^^'^'^ involves a conversion of the cc pair into 
the rj' via two gluon exchanges. Although the charm content of the 1]' is a priori expected to 
be small, its contribution is potentially important because the CKM mixing angle VcbV*^ is 
of the same order of magnitude as that of the penguin amplitude [cf. Eq. (5.6)] and yet its 
effective coefficient 02 is larger than the penguin coefficients by an order of magnitude. The 
decay constant /^,, defined by (0|c7^75c|?7') = i/^/?/*, has been estimated to be = (50—180) 
MeV, based on the OPE, large- A^c approach and QCD low energy theorems [39]. It was 
claimed in [39] that \f^, \ ~ 140 MeV is needed in order to exhaust the CLEO observation of 
B"^ rj'K^ and B rj' + X hy the mechanism 6 cc + s 77' + s via gluon exchanges. 
However, a large value of f^, seems to be ruled out for several reasons [50]. For example. 



§§A two-mixing-angle parametrization of the rj and r]' wave functions: r]' = ris sin 6*8 + rjo cos Oq, rj = 
rjs cos 9s — rjo sin ^o, is employed in [24] for the calculation of S — > r)'{r])K. However, in the absence 
of mixing with other pseudoscalar mesons, this parametrization will destroy the orthogonality of 
the physical states rj and rj' if Oq ^ 6s- Due to SU(3) breaking the matrix elements (0|A^^^^ 1 778(0)) 
do not vanish in general and they will induce a two-angle mixing among the decay constants: 

Based on the ansatz that the decay constants in the quark flavor basis follow the pattern of particle 
state mixing, relations between 9s, 9q and 9 are derived in [53], where 9 is the rj — rj' mixing 
angle introduced in (5.11). It is found in [53] that phenomenologically 9s = —21.2°, = —9.2° 
and 9 = —15.4°. It must be accentuated that the two-mixing angle formalism proposed in [54,53] 
applies to the decay constants of the rj' and rj rather than to their wave functions. Numerically, we 
find that the branching ratios shown in Table I (see below) calculated in one-angle and two-angle 
mixing schemes are different by at most 7%. In the present paper we shall employ the former 
scheme. 
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from the data of J/ip rjc'y and J/ip ^ rj'j, one can show that |/^^, | > 6 MeV, where 
the lower bound corresponds to the nonrelativistic quark model estimate. Based on the 777 
and T)'^ transition form factor data, the range of allowed f^, was recently estimated to be 
—65 MeV < /^/ < 15 MeV [55]. A most recent reevaluation of f^, along the hne of [39] yields 
= —(12.3 ~ 18.4) MeV [56], which is in strong contradiction in magnitude and sign to 
the estimate of [39]. The sign of f^, can be fixed by using QCD anomaly and is found to 
be negative [47] (see also [49,53,56]). In the presence of the charm content in the 770, an 
additional mixing angle 9c is needed to be introduced: 

I770) = — F cos 9c\uu + dd + ss) + sin 9c\cc), 



1 

73' 



|?7c) = 1= sm9c\uu + dd + ss) + cos 9c\cc). (5-17) 



Then f^^, = cos 9 tan 9c fr,^ and = — sin 9 tan ^c/r?c > where the decay constant /^^ can be 
extracted from rjc — > 77, and 9c from J/ip ^ rjcj and J/ip ^ rj'j [24]. In the present paper 
we shall use 

f^, = -6 MeV, = - tan 9f^, = -2.4 MeV, (5.18) 

for 9 — —22° (see below) , which are very close to the values 

f^, = -(6.3 ± 0.6) MeV, = -(2.4 ± 0.2) MeV (5.19) 

obtained in [53]. 

Using F(f^(0) = 0.34 [32], y^Fo'^yO) = 0.33 *** for form factors, fo = = 
ms(l GeV) = 150 MeV, 9 = —19.5° and Eq. (5.10) for the matrix element {ti'\s^5s\0), we find 
that B{B — > rj'K) = (0.9—1.0) x 10^^ and it is insensitive to A'^^*' and the choice of Wolfenstein 
parameters p and rj so long as Vp^ + v'^ ~ 0.37, where Nf = Nf{V+A) = Nf{V-A). The 
discrepancy between theory and experiment can be greatly improved by the accumulation of 
several enhancements. First of all, the running quark masses appearing in the {S — P){S+P) 
matrix elements should be applied at the scale fi = (9(mf,) as given in Eq. (2.13) so that the 
{S — P){S + P) matrix element is enhanced due to the decrease of ms{fi) at fi = nii,. (The 
sensitivity of the branching ratio to was first noticed in [42].) Second, a recent analysis 
of the data of r], rj' — >• 77 and r], rj' — >• ttttj yields [57] 

y- = 1.38 ±0.22, ^ = 1.06 ±0.03, ^ = -22.0° ± 3.3°, (5.20) 

Jit Jit 



***The form factors Fff '(0) = 0.254 and Fq '\0) = 0.307 given in [7] do not take into account 
the wave function normalization of the physical r]' and tj states. Since it is not clear to us what is 
the r) — r)' mixing angle employed in [7], we shall follow [24,27] to use the nonet symmetry relation 
V3F(f'?o(0) = VeF^^^iO) = ^(^^^(0) ^ 0.33 to obtain F(f^°, F(f''« and hence the form factors F^''' 
as well as Fq^ for a given 6. 
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implying some SU(3) breaking in the decay constants. Applying the new values of the 
aforementioned parameters, the result for the branching ratio of — >■ r]'K'^ is shown in 
Fig. 6 vs l/N^^ (see the lower set of solid and dotted curves). We find that B{B^ — > r]'K^) 
is enhanced from (0.9 — 1.0) x 10~^ to (2 — 3) x 10~^. The latter result is in agreement with 
[24] (see the lower set of curves with negative f^, in Fig. 17 of [24]). The enhancement is due 
mainly to the running strange quark mass at ^ = mb and SU(3) breaking effects in the decay 
constants /q and /§. From Fig. 6 we see that (i) in the absence of the anomaly contribution 
to {ti'\s^5s\0), the branching ratios (the upper set of solid and dotted curves) will be further 
enhanced in a sizable way (of course, it is erroneous to neglect such an anomaly effect), and 
(ii) the contribution of cc conversion into the f]' becomes destructive when < 0.28. 

This is understandable because 02 becomes negative at small values of l/N^^ so that the 
term 02^^^^'^'-' contributes in opposite sign to the penguin amphtudes. Therefore, the charm 
content of the ri' is not welcome for explaining B{B — > rj'K) at small l/Nf^. 

Thus far it has been assumed in the analysis of B ^ rj'K that the nonfactorizablc effects 
lumped into Oj via (iV^^)j are the same for z = 1, 2, ■ • • , 10. However, we have pointed out in 
Sec. Ill that N^^ {V — A) in hadronic charmless B decays is most likely very similar to that in 
B Dn, namely N^^iV — A) ~ N^^{B Dn) 2. In fact, we just showed that the charm 
content of the rj' will make the discrepancy between theory and experiment even worse at 
small values of 1/Nf if Nf{V - A) is the same as Nf{V + A). Setting Nf{V - A) = 2, 
we find that (see Figs. 7 and 8) the decay rates of 5 — > rj'K are considerably enhanced 
especially at small l/Nf{V + A). That is, B{B^ rj'K^) at l/Nf{V + A) < 0.2 is 
enhanced from (2.5 — 3) x 10~^ to (3.7 — 5) x 10'"''. First, the 1]' charm content contribution 
a2-^c^^'''''* now always contributes in the right direction to the decay rate irrespective of the 
value of N^^{V + A). Second, the interference in the spectator amplitudes of B^ — > rj'K^ is 
constructive. Third, the term proportional to 

2(a3 - a5)Xf + {a, + a, - 05)^^ ^'^'^ (5.21) 

in Eqs. (5.3) and (5.4) is enhanced when {Nf)^ = {Nf)4 = 2. It is evident from Fig. 8 
that the measurement of B^ — > rj'K^ is well explained in the present framework based on 
the Standard Model within the allowed range l/Nf{V + A)< 0.23 extracted from B^ 
(j)K^. Contrary to some early claims, we see that it is not necessary to invoke some new 
mechanisms, say the SU(3)-singlet contribution S' [43], to explain the data. The agreement 
with experiment provides another strong support for N^^[V — A) ~ 2 and for the relation 
Nf{V + A)> Nf{V - A). As for the decay B^ rj'K^, the predicted branching ratio, 
say 4 X 10~^ at our preferred value Nf^{V + ^4) ~ 5 (see Table I), is compatible with the 
data, though it is on the lower side. For a slightly enhanced decay constant Ki —15 MeV, 
as implied by a recent theoretical estimate [56], we obtain B{B — > rj'K) = (4.6 — 5.9) x 10~^ 
at l/Nf{V + A)< 0.2, which agrees with experiment very nicely. Note that the CLEO 
data of B^ —> rj'K^ and B^ rj'K are in good agreement within one sigma error [see (5.2)], 
though the charged mode is more reliable. It is conceivable that when errors are improved 
and refined, the two values will converge eventually. 

We have also studied the decays B — > r]K, rj'K* , rjK* . The decay amplitude of S — > riK 
is the same as S — > rj'K except for a trivial replacement of the index 77' by 77. As a general 
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rule, the factorizable amplitude of -B ^ rj^'^K* can be obtained from the B rj^'^ K one 
by (i) replacing the term mp/[(mi + m2)(m3 — rrii)] by — mp/[(mi + m2){mz + 777.4)] and 
the index K hy K*, and (ii) discarding the [S — P){S + P) contribution associated with 
Pqj. example, the decay amplitude of B" — > r]'K*~ can be easily read from (5.3) 

to be: 



It/ 111 

(a4 + aio)X^^'''''^'^ + U3 + a4 - 05 + -07 - -ag - -aio 



+ (203 - 2a5 - ^07 + ^ag) Xf + (03 - as - a7 + ag)^^ 



+ ( a4 + aio - 2(a6 + a^) , ^ "f' , ) X^^'^'^*^ 

\ [rUs + 777„) (777b + "^u) / 



(5.22) 



with 



XiBr,',K*) ^ (X*-|(,-^)^_JO)(77'|(tZ6),_JS-) 

= -2fK*mK*Fi'' {mj^*){e ■ pb), 
XiBK*,,') ^ (^'l^qq)^_Jo)^K*-\i7sb),_jB-) 

j^iB,r,'Kn ^ (^r,'K*-\{-su),_,\0){0\{ub),_jB-). (5.23) 

From Table I we see that the electroweak penguin is generally small due to the smallness 
of its Wilson coefficients, but it does play an essential role in the decays B^ rjK^ and 
B^ — > TjK'^. It is interesting to note that the branching ratios of i? — > 7]^'^K'^*^ are all 
less than 1 x 10~^ except for B r)'K, which has a very large branching ratio, of order 
(4 — 6) X 10~^. It has been argued in [39] that B{B — > rj'K*) is about twice larger than that 
of S — > rj'K, which is certainly not the case in our calculation. The ratios of various decay 
rates are predicted to be 

B{B rj'K) _ r 72 B{B rj'K*) _ ( 0.06 charged B; 

B{B r)K) ~ I 296 ' B{B r)K*) ~ I 0.02 neutral S, ^ ' ' 

for positive p. The destructive (constructive) interference between the terms 
a^X(BK,r,^ ') expl ains the ratio B{B — > ri'K)/B{B — >• rjK): X^^^'^'^ has a sign opposite to 
X^BK,ri) Qjjg easily see from the wave functions of the 77 and 77', Eq. (5.11). Since 
the sign of a^X'^^^* ''^'^ is flipped in S ^ rj'^'^K* decays, the interference effect becomes the 
other way around: constructive in S — > rjK* and destructive in S — > rj'K*. 

To discuss the decays B — > 7;*^'^7r(p), we consider B" — > rj'n' as an illustration. Its decay 
amplitude is 
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Table I. Branching ratios averaged over CP-conjugate modes for charmless B decays to the rj' and rj, where 
"Tree" refers to branching ratios from tree diagrams only, "Tree+QCD" from tree and QCD penguin dia- 
grams, and "Full" denotes full contributions from tree, QCD and electroweak (EW) penguin diagrams in 
conjunction with contributions from the process cc — > jjq- Predictions are for fc^ = r] = 0.35, p = —0.12 

(the first number in parentheses) and r] = 0.34, p = 0.16 (the second number in parentheses). The decay 
constants f^, = —6 MeV and = —2.4 MeV are used. The effective number of colors is taken to be 
N^^{V — A) =2 and N^^{V + A) = 5. The running quark masses at /i = m6 are given by (2.13). 
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where 
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and 



(5.26) 



^ (V|(g-g),_jO)(7r-|(J6),_jS-) = i/^M^B - ml)Ft{ml,), 
XiB,.'^) ^ {r,'n-\{du),_jO){0\{vb),_,\B-). (5.27) 

In deriving (5.25) we have applied the matrix elements 

{V'\ui5u\0) = {v'\dj5d\0) = V W\si5s\0), (5.28) 

with r^/ being given by (5.26). 

Since VubV*^, VcbV*^, VtbV^^ arc all comparable in magnitude [cf. Eq. (3.15)] and since the 
Wilson coefficients of penguin operators are rather small, it is expected that B — r^'^'V, r]'^'^ p 
are dominated by spectator diagrams^^^. Prom Table I we see that this is indeed the case 
except for the decay modes B^ — > jy^'^Tr^ which are penguin dominated. To compute the 
decay rate of S — > r]T:{p) we have applied the matrix element {r]\u-)^u\Qi) — r^(77|s75s|0) with 



_ ^ y 2/o^ - /I cosg - v/2sin6' ^c: 9Q^ 

The mechanism of cc — > 770 is less significant m. B ^ ri^'^7r{p) decays because it does not 
gain advantage from the quark mixing angle as in the case of S — > rj^">K{K*). We see from 
Table I the minor role played by the charm content of the 77' except for the decay S° — > rj'n^. 
In general, the decay rates of S — > r]^'^7r{p) are not sensitive to the values of Nf^{V — A) 
and Nf^(V + A) and do not vary significantly from channel to channel: 

B[B^ 7]^'^7r{p)] - (3 - 10) X 10-^ B[B^ V^'^^^ip)] (0-2 - 4) x 10"^. (5.30) 

It is interesting to note that B[B — > 777r(p)] > B[B — > ri'7r{p)]. 



^^^I^Thc matrix element {r]'\wy5u\0) can be obtained from [53] and it is slightly different from the 
corresponding one in [52,51]: 

/8 cos + 4j/o sin 6* i 
/8Cos6l- V2/osm6' 2/^ 



M^The branching ratios of S ^ rj^'^ir, rj^'^p are largely overestimated in [40,50] as the incorrect 
matrix element {rj' {u'y^ulO) = —im'^,f^j/{2mu) is applied there. 
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VI. DIFFICULTIES WITH B' R-oj 

Up to now we have shown that CLEO results on hadronic charmless B decays can be 
satisfactorily explained provided that Nf{V - A) ^ 2 and Nf{V + A) ^ 0{A). However, 
there is one CLEO measurement, namely the decay — > cu/T^, that is beyond our expla- 
nation and hence may impose a potentially serious difficulty. In this Section we will first 
explore the problem and then proceed to suggest some possible solutions. 

The decay amplitude of B~ —>■ ujK~ is very similar to B~ u!Tr~ and has the expression 



ts 



04 + aio — 2(a6 



{rrib + niu) {rUs + m^) ) 



+ ^(4a3 + 4a5 + a7 + a9)Xf^'") 

\ (mfe + niu) [ms + my) J 



(6.1) 



where 



^'-) = {u\{uu),_,mK-\{sh),_,\B-) = -v^^m^Ff^(m^)(£ -ps), 

^ (a;X-|(.-^),_jO)(0|(fZ6),_j5-). (6.2) 

We see from Fig. 9 that the calculated branching ratio using Nf^ (V — A) — 2, Ff^{Q) — 0.34 
and Aq'^{0) = 0.28/\/2 [7] is too small compared to experiment [34]: 

B{B^ uK^) = (l.5t°:^ ± 0.2) X 10-^ (6.3) 

In fact, all the region of l/Nf{V + A) < 0.9 is excluded. Nevertheless, if Nf{V - A) 
is taken to be the same as Nf^(y + A), then a rather small value of '\-/N'^^ < 0.05 is 
experimentally allowed [24,27] (see Fig. 10). In other words, N^^ is preferred to be very 
large in B^ —>■ ujK^. In our opinion, however, a very large value of Nf^iV — A) is rather 
unhkely for several reasons: (i) A small Nl^{V — A) pa 2 is favored in other charmless 
B decays: B — > tttt, tto; and B — > rf'K. (ii) It will lead to a too large nonfactorizable 
term, which is not consistent with the small nonfactorizable effect observed in the spectator 
amplitudes ol B ^ Dtt and the picture that the nonperturbative feature of nonfactorizable 
effects is loose in the energetic two-body decays of the B meson, as we have elaborated 
before (see the end of Sec. IV). It thus appears to us that the observed large decay rates of 
B^ — > ljK^ is attributed to other mechanisms rather than to a very large value of N^^. 

So far we have neglected three effects in the consideration of B"^ — > cuK^: VT-annihilation, 
space-like penguin diagrams and final-state interactions (FSI); all of them are difficult to 
estimate. In order to understand why B{B'^ — > lott"^) <^ B{B^ — > uK"^) experimentally, we 
need a mechanism which will only enhance the latter. It appears that FSI may play this 
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role. Since — > ujK~ involves only a single isospin amplitude, inelastic scattering will 
be the dominant effect of FSI. For example, b — >■ ccs and b — > uus modes can mix with 
each other so that the decay B~ — > ujK" arises either from b — > ccs or indirectly through 
B- D^D*- or D*°D- (via b ccs) with a rescattering D^D*- (or D*^D~) uK~. For 
the decay B" — > con', the inelastic scattering B~ — > {DD*} — > con' is Cabibbo suppressed. 
Therefore, it is possible that B~ ujK~ receives large FSI from inelastic scattering but 
B~ UJTT~ does not. Since B^ uK^ does not receive contributions from 14^-annihilation, 
its measurement can be used to test the relative strength between FSI and annihilation 
terms. If the branching ratios of B^ — > ujK^ and B^ — > ujK^ are close, this will imply the 
importance of FSI. 

VII. DISCUSSIONS AND CONCLUSIONS 

For a given effective weak Hamiltonian, there are two important issues in the study of the 
hadronic matrix elements for nonleptonic decays of heavy mesons: one is the renormaliza- 
tion scale and scheme dependence of the matrix element, and the other is the nonfactorizable 
effect. For the former, we have emphasized that it is important to first evaluate the vertex 
and penguin corrections to the matrix element of 4-quark operators at the scale so that 
{0{fi)) = g{^){0) tree and then apply factorization or any model calculation to (O)trcc- The 
resulting effective coefficients cf^ = Ci{fi)g{fi) are renormalization-scale and -scheme indepen- 
dent. We pointed out that while c^|| ~ (^i^if^) at /x = mb{mh) for current-current operators, 
the real parts of c|^g are about one and half times larger than the leading-order penguin 
Wilson coefficients. This means that to describe the hadronic charmless B decays domi- 
nated by penguin diagrams, it is necessary and inevitable to take into account the penguin 
corrections to the 4-quark operators. 

Nonfactorizable effects in hadronic matrix elements of B ^ PP, VP decays can be 
parameterized in terms of the effective number of colors N^^ in the so-called generalized 
factorization scheme; the deviation of l/Nf^ from 1/Nc {Nc = 3) characterizes the nonfac- 
torizable effect. We show that, contrary to the common assumption, N^^{V + A) induced 
by the {V — A){V + A) operators 05,6,7,8 are theoretically and experimentally different from 
Nf{V - A) generated by the {V - A){V - A) operators. The CLEO data of B^ un^ 
available last year clearly indicate that Nf^{V — A) is favored to be small, Nf^{V — A) < 2.9 . 
This is consistent with the observation that N^^(y — /I) pa 2 in i? ^ Dn decays. Unfortu- 
nately, the significance of B^ ujtt^ is reduced in the recent CLEO analysis and only an 
upper limit is quoted. Therefore, a measurement of its branching ratio is urgently needed. 
In analogue to the class-Ill B Dn decays, the interference effect of spectator amplitudes 
in charged B decays B~ — > tt^tt^, p"7r°, n^p^ is sensitive to N^^(V — A); measurements of 
them [see (3.23)] will be very useful to pin down the value of N^^{V — A). 

Contrary to the nonfactorizable effects in spectator-dominated rare B decays, we found 
that N^^{V + A) extracted from the penguin-dominated decay — > (pK^^ is larger than 
N^^iy — A). This means that nonfactorizable effects in tree and penguin amplitudes behave 
differently. It turns out this observation is the key element for understanding the CLEO 
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measurement of B rj'K. In the conventional way of treating A^^^(V" + A) and Nf^iV — A) 
in the same manner, the branching ratio of — >■ r]'K'^ after including the anomaly effect 
in the matrix element (77'|s75s|0) is naively only of order 1 x 10""^. The running strange 
quark mass at = and SU(3) breaking in the decay constants /s and /o will enhance 
B{B rj'K) to the order of (2 - 3) x IQ-^ with f^, = -6 MeV. This is still lower than the 
central value of the CLEO measurements. Also, the charm content of the rj' is not welcome 
for explaining the decay rate of -B ^ rj'K at small values of We showed that the 

fact that N^^{V + A) > N^^{V — A) ^ 2 will substantially enhance the branching ratio of 
B^ ri'K^ to (3.7-5) x 10^^ at 1 /NfiV+A) < 0.2. Unlike the previous analysis, the small 
charm content of the rj' is now always in the right direction for enhancement irrespective 
of the values of l/Nf{V + A). The predicted branching ratio of 5° rj'K^ IS m good 
agreement with experiment and the calculation of B"^ — > rj'K^ is compatible with the data. 
For a slightly enhanced —15 MeV, as implied by a recent theoretical estimate, we 

found that the agreement of the predicted branching ratio for B ri'K with experiment is 
very impressive. It is thus important to pin down the decay constant /^,, recalling that the 
commonly used value |/^^,| = 6 MeV is extracted from experiment within the nonrelativistic 
quark model framework. We conclude that no new mechanism in the Standard Model or 
new physics beyond the Standard Model is needed to explain B rj'K. We have also 
analyzed charmless B decays into the rj' and 77 in some detail. The branching ratios of 
the spectator-dominated decays B — > r]^'^n, r]^'^p were largely overestimated in the previous 
analysis because the matrix element (rj^'^u'j^ulO) was not evaluated correctly before. 

Although the CLEO measurements of hadronic charmless B decays are satisfactorily 
explained in the present framework, we found that it is difficult to understand the experi- 
mental observation that r{B^ un'^) < r{B^ uoK^). The calculated branching ratio of 
B^ — s> uoK"^ is too small compared to experiment. We conjecture that final-state interactions 
via inelastic scattering may contribute in a sizable way to B^ uK^, but are negligible for 
B^ — > cuTT^ due to the Cabibbo-angle suppression. Clearly this decay mode deserves further 
serious investigation and a measurement of the neutral decay mode S° — > cuK'^ will be very 
useful to clarify the issue. 

Under the factorization hypothesis, the decays B (j)K and B (l)K* should have 
almost the same branching ratios, a prediction not borne out by current data. Therefore, it 
is crucial to measure the charged and neutral decay modes of B 4'{K, K*) in order to see 
if the generalized factorization approach is applicbale to B ^ (pK* decay. 

To conclude, based on the available CLEO data on hadronic charmless two-body decays 
of the B meson, we have shown that the nonfactorizable effect induced by the (V — A){V + A) 
operators is different from that generated by the {V — A){V — A) operators. This is the key 
element for explaining the CLEO measurement of S — > ij'K. 
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Figure Captions 



Fig. 1. The branching ratio of B"^ — > cvn"^ vs l/Nf^. The sohd and dashed curves are 
for T] = 0.34, p — 0.16 and rj — 0.35, p — —0.12, respectively. The solid thick hnes are the 
CLEO measurements with one sigma errors. 

Fig. 2. Same as Fig. 1 except that the branching ratio is plotted against l/Nf^{V + A) 
with Nf{V - A) being fixed at the value of 2. 

Fig. 3. Same as Fig. 1 except for — > 7r''=7r°. The thick dotted line is the CLEO upper 
limit [see (3.20)]. 

Fig. 4. The branching ratio of (pK^ vs for r] = 0.34 and p = 0.16. The 

dotted curve is for N^^(y + A) = N^^{V — A) = N^^ and the solid curve is the branching 
ratio against l/Nf{V + A) with Nf{V - A) being fixed to be 2. The solid thick hue is the 
CLEO upper limit. 

Fig. 5. Same as Fig. 1 except for B (pK*. 

Fig. 6. The branching ratio of B^ rj'K^ as a function of for 7] = 0.34 and 

p = 0.16. The charm content of the 1]' with f^, = — 6 MeV contributes to the sohd curves, 
but not to the dotted curves. The lower set of sohd and dotted curves takes into account the 
anomaly contribution to (?7'|s75s|0) [see Eq. (5.10)], whereas the upper set does not. The 
solid thick lines are the CLEO measurements with one sigma errors. 

Fig. 7. Same as Fig. 6 except that the branching ratio is plotted against 1/N^^{V + A) 
with N^^iy — A) being fixed at the value of 2. The anomaly contribution to (?7'|s75s|0) is 
included. 

Fig. 8. Same as Fig. 7 except for S° rj'K'^. 

Fig. 9. The branching ratio of B^ uK^ vs l/Nf{V + A) with Nf{V - A) being fixed 
to be 2. The sohd and dashed curves are for 77 = 0.34, p — 0.16 and 77 = 0.35, p — —0.12, 
respectively. The sohd thick hnes are the CLEO measurements with one sigma errors. 

Fig. 10. Same as Fig. 9 except that Nf{V + A) = Nf{V - A) = Nf. 
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